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INTRODUCTION 

1 2 

The  NASA  “PAC  2"  program  for  calculation  of  thermodynamic 
quantities  of  chemical  species  has  as  one  of  its  options  the  calculation 
of  heat  capacity,  enthalpy,  and  entropy,  and  the  fitting  of  these  to  a 
polynomial  function  of  the  temperature.  The  polynomial  coefficients  may 
be  specified  to  be  in  the  form  suitable  as  input  data  for  the  companion 
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NASA  program  to  calculate  thermodynamic  equilibrium.  This  latter 
program  has  found  considerable  use  at  NAVWPNSUPPCEN  Crane. 

In  expanding  the  number  of  elements  which  may  be  involved  in  the 
equilibrium  program  It  has  been  necessary  to  calculate  the  polynomial 
coefficients  for  many  chemical  species,  using  PAC  2. 

In  order  to  get  accurate  polynomial  coefficients  it  is  necessary  to 
furnish  thermodynamic  data  at  many  closely  spaced  temperatures  as  input 
to  the  PAC  2 program.  This  Involves  considerable  keypunching.  However 
suitable  data  for  a wide  variety  of  chemical  species  are  in  the  JANAF 
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tables.  Therefore,  we  have  modified  PAC  2 to  read  the  data  from  the 
JANAF  tape  rather  than  from  cards. 

Another  time-saving  modification  was  related  to  the  requirement  of 
the  PAC  2 program  for  input  data  at  the  transition  points.  This  is  not 
available  In  the  JANAF  tables,  and  it  had  been  necessary  to  obtain  it  by 
hand  calculation.  The  modified  program  now  does  this  by  extrapolation 
from  the  other  data  furnished. 

GENERAL  DESCRIPTION  OF  PROGRAM 

The  PAC  2 program  has  options  for  a number  of  different  types  of 
functions  to  be  calculated.  In  a particular  run  a list  of  species  may  be 
given,  with  a different  option  requested  for  each.  For  any  particular 
species  a variety  of  methods  are  available,  each  requiring  a different 
type  of  Input  data. 

Typically,  using  the  JANAF  tables**  as  a data  source,  one  would  read 
in  the  thermodynamic  functions  at  100  degree  intervals  for  the  solid  and 
some  of  the  liquid  phases  using  the  READIN  method.  For  liquid  phases 
where  the  heat  capacity  is  giver  as  constant,  or  can  readily  be  seen  to 
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be  a linear  function  of  temperature,  it  is  most  efficient  to  furnish  this 
information  and  use  the  COEF  method.  For  gaseous  species  one  of  the 
theoretical  methods  (JANAF,  PANDK,  etc.)  using  spectroscopic  data 
required  the  least  effort  in  preparing  input. 

The  new  option  was  designed  to  take  the  thermodynamic  data  from  the 
magnetic  tape  using  the  READIN  method,  and  then  use  the  LSTSQS  option  to 
obtain  the  polynomial  coefficients.  It  is  probably  possible  to  change 
options  and  methods  from  one  species  to  the  next  during  a given  run  since 
the  other  options  and  methods  were  not  changed.  This  was  specifically 
demonstrated  using  the  JANAF  method  to  generate  the  polynomial 
temperature  coefficients  of  the  thermodynamic  functions  for  one  species 
in  a list,  while  using  READIN  from  tape  for  the  others.  It  is  also 
possible  to  turn  off  the  tape  and  furnish  the  input  for  the  READIN  method 
from  cards  for  arbitrary  species  within  a list.  It  is  not  possible  to 
switch  between  tape  and  cards  as  input  for  the  READIN  method  for 
different  phases  of  a particular  species.  It  is  also  not  possible  to 
change  methods  on  going  from  one  phase  to  the  next  with  the  tape  on 
option,  for  Instance  to  use  READIN  for  the  solid  phases,  then  use  COEF 
for  the  liquid  phase.  However,  this  change  of  method  can  still  be  done 
If  READIN  is  performed  using  data  from  cards,  with  TAPEOFF  specified. 

Despite  this  latter  restriction,  the  simplification  of  the 
preparation  of  input  for  reading  data  from  tape  more  than  compensates  for 
the  extra  computer  time  involved  in  doing  least  squares  fits  to  phases 
with  constant  heat  capacities. 

SPECIFIC  INPUT  CARDS 

The  cards  to  be  used  for  each  species,  with  changes  made  from  the  PAC 
2 program  as  described  in  references  1 and  2,  are  as  follows: 

1.  Formula  Card  - no  change  from  reference  1. 

2.  TAPEOFF  Card  - requires  the  data  to  be  read  from  cards. 

In  the  absence  of  this  card  data  are  read  from  the  tape. 

3.  LSTSQS  - no  change  from  reference  1. 

4.  Method  Card  - contains  melting  point  or  transition 
temperature  and  heat  of  transition,  as  in  reference  1.  The  only  change 
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Is  that  at  a melting  point,  but  not  at  other  phase  transitions,  the  new 
heat  of  formation  must  be  provided,  since  subsequent  data  from  the  tape 
are  relative  to  It. 

5.  Data  Card  - a single  card  gives  the  species  number  (see  a 
listing  of  the  JANAF  data  tape)  and  the  temperature  range,  in  even 
hundreds,  for  the  new  phase,  in  314  format.  The  method  and  data  cards 
are  repeated  for  each  phase,  as  in  reference  1. 

6.  Finish  Card  - as  in  reference  1. 

OPERATION  OF  THE  MODIFIED  PROGRAM 

Upon  reading  the  data  card  the  program  searches  the  data  tape  for  the 
species  number  given,  and  there  reads  Cp,  S,  and  H-H^gg  data  at  each 
100  degree  interval  between  the  limits  given  on  the  data  card.  The  next 
Input  card  Is  then  read,  giving  the  next  transition  temperature.  The 
program  then  computes  the  polynomial  coefficients  by  the  least  squares 
procedure;  and  these  coefficients  are  then  used  to  calculate  the 
thermodynamic  functions  Cp,  S,  and  H-H^gg  for  the  low  temperature 
phase  at  the  transition  temperature.  The  enthalpy  of  transition  is  then 
used  to  calculate  the  enthalpy  and  entropy  for  the  high  temperature  phase 
at  this  transition. 

The  next  data  input  card  is  read,  and  the  program  then  selects  from 
the  data  tape  the  values  of  Cp,  S,  and  H^gg  above  the  transition. 

In  the  case  the  transition  is  a melting  point,  the  new  H-H^gg  data  are 
corrected  to  refer  to  the  heat  of  formation  of  the  solid  at  298.15°. 

The  heat  capacity  of  the  high  temperature  phase  at  the  transition  is  then 
obtained  by  backward  extrapolation  from  the  succeeding  1 to  3 (as 
available)  values  of  Cp  . The  program  now  has  all  the  data  needed  for 
the  least  squares  procedure  to  evaluate  the  polynomial  coefficients  for 
this  phase.  The  procedure  is  repeated  until  all  requested  phases  have 
been  calculated. 

The  input  information  from  cards  and  tapes,  and  the  smoothed  values, 
are  printed.  The  polynomial  coefficients  are  put  onto  a punch  tape  in 
the  format  required  for  the  data  input  cards  of  reference  3.  This  is 
just  as  is  done  In  reference  1.  The  data  calculated  for  the  transition 
temperature  are  Included  in  the  printout. 
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APPENDIX  A 


ASDH 

CHS 

DD 

ICODE 

IFIRST 

IJW 

IKODE 

IOLDC 

I TAPE 

NFRST 

NR  EC 

PHASDW 

PHAST 

POLTU 

POLTUL 

TAPE8 

TEST  (2) 

TLAGR 

TLL 

TUL 

TT 

TTON 
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New  Variables 

Correction  term  to  ASINDH  for  new  phase 
Contents  of  a card 

Array  of  thermodynamic  functions  from  tape 

Code  number  of  current  species 

Next  tape  access  will  be  the  first  for  this  species 

Buffer  area 

Code  number  on  card  from  tape 

Code  number  of  last  species  obtained  from  tape 

Number  of  temperatures  selected  from  tape  for  input  to  the 

least  squares  routine 

Flag  to  indicate  that  this  is  the  first  card  for  the  current 
species 

Number  of  temperatures  read  from  tape 
Increment  in  enthalpy  at  phase  transition 
Phase  transition  temperature 
Lower  temperature  limit  for  previous  polynomial 
Upper  temperature  limit  for  previous  polynomial 
Logical  file  name  for  thermodynamic  data  tape 
Input  thermodynamic  data  from  tape  (if  true) 

Perform  a Lagrangian  interpolation  (if  true) 

Lower  temperature  limit  for  data  to  be  selected 
Upper  temperature  limit  for  data  to  be  selected 
Array  of  temperatures  from  tape 
Thermodynamic  data  is  to  be  read  from  tape  (if  true) 
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